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Abstract. It is known that any subharmonic quadrature domain in two di- 
mensions satisfies a natural inner ball condition, in other words there is a 
specific upper bound on the curvature of the boundary. This result directly 
applies to free boundaries appearing in obstacle type problems and in Hele- 
Shaw flow. In the present paper we make partial progress on the corresponding 
question in higher dimensions. Specifically, we prove the equivalence between 
several different ways to formulate the inner ball condition, and we compute 
the Brouwer degree for a geometrically important mapping related to the 
Schwarz potential of the boundary. The latter gives in particular a new proof 
in the two dimensional case. 
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1. Introduction 

In the present paper we study the curvature of a free boundary which comes up in 
some obstacle type problems, more specifically Laplacian growth (Hele-Shaw flow 
moving boundary problem), quadrature domains and partial balayage. The final 
aim of the investigations is to show that the free boundary in question satisfies a 
certain inner ball condition (a specific upper bound on the curvature). This goal 
has previously been achieved in the case of two dimensions (see [7] and [5] ; compare 
also [H]). Here we shall give some partial results (but no complete solution) in 
higher dimensions, and in passing also obtain a new proof for the two dimensional 
case. 
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The geometric property we aim at proving can most easily be stated in terms 
of quadrature domains for subharmonic functions jlOj . Let /u, a positive Borel mea- 
sure with compact support in R™. A bounded open set f2 C R n is called a quad- 
rature open set for subharmonic functions with respect to fj, if fx(JS, n \ f2) = 
and 

/ hd/i < / h dm 
Jn Jn 

for all integrable (with respect to Lebesgue measure dm) subharmonic functions h 
in tt. A quadrature open set which is connected is a quadrature domain. There is a 
natural process of balayage of measures to a prescribed density (partial balayage) 
by which fi can be constructed from fi when it exists, see e.g. [TU], [BJ, [5]. This 
balayage process is also equivalent to solving a certain obstacle problem [11] . 

In terms of the difference u = — U n between the Newtonian potentials of 
H and f2, the latter considered as a body of density one, the quadrature property 
spells out to 

'u>0 inR™, 
< u = outside Q, 
Au = X n-Li hiK". 
The function u appearing here is sometimes called the modified Schwarz potential 
of dfl, see [15] for example. What has been proved in two dimensions, and what we 
like to extend to higher dimension, is that fl in the above situation can be written 
as the union of open balls centered in the closed convex hull K of supp /i: 

n = U xeK B(x,r{x)). (1.1) 

Here r(x) > denotes the radius of the ball at x, allowing the possibility r(x) = 0, 
i.e., that the ball is empty. We refer to (jl.ll) as ft satisfying the inner ball property 
with respect to K. For further discussion and motivations in the present context, 
see H, 0. 

The inner ball property concerns the geometry of dfl outside any closed half- 
space H containing supp [i. By a certain "localization" procedure the part Sl\H 
of Q which is outside H can be shown to be identical with a quadrature open set 
for some positive measure with support on dH, see [6J, [13] . For this reason it is 
enough to prove (jl .![) in the case that jj, has support in a hyperplane, and by a 
further localization one may even assume (i to have a continuous density on it. For 
convenience we shall take the hyperplane in question to be {x € R™ : x n = 0}. It is 
known (see [10] , [6J ) that given any positive measure with compact support in this 
hyperplane, which we identify with R™ -1 , there is a uniquely determined quadra- 
ture open set, which moreover is symmetric about the hyperplane and is convex in 
the a;„-direction (i.e., the intersection with any straight line perpendicular to the 
hyperplane is connected). Thus the quadrature open set can be described in terms 
of a graph of a function g defined in an open subset D of the hyperplane. It is 
known from the regularity theory of free boundaries (see [4] , pQ , and in the present 
context [6j) that this function g is real analytic. It will be enough to consider the 
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case that D is connected, because in the disconnected case the discussions will 
apply to each component separately. 

In the paper we shall therefore study the geometry of quadrature domains 
for a measure with support in the hyperplane R™ _1 , the corresponding modified 
Schwarz potential u, and various differential geometric objects derived from it. 
The paper consists of two main parts. The first part starts with some differential 
geometric preliminaries (Section [2J and ends up with a proof of equivalence of 
several different formulations of the inner ball condition (Section [3]). This part is 
analogous to a corresponding part in [7], but new difficulties appear in the higher 
dimensional case. In the second part of the paper, Section 21 we develop tools for 
studying the geometry by means of vector fields and differential forms defined in 
terms of u. 

To briefly explain, in terms the two dimensional situation, what we do in the 
second part of the paper, let fi + denote the part of the quadrature domain which 
lies in the upper half space (half plane) and let S(z) be the Schwarz function (see 
[2J, [TS]) of (dfl) + . This is in our case given by S(z) = z — 4^, it is analytic in 
Q + , and equals z on (9f2) + . What we study in the second part of the paper is the 
higher dimensional counterpart of the mapping a : Q + — > C defined by 

1 1 2 du du 

a{z) = -zS{z) = -r -r^+i^, 

the last member referring to polar coordinates. We also study in higher dimensions 
that curve 7 which in two dimensions is defined by |p = 0. In the two dimensional 
case a proof of the inner ball property can be based on the topological property of 
7 that it separates the two domains < and > from each other, see [8]. 
There seems to be no direct counterpart of this kind of proof in higher dimensions. 

A slightly different proof in two dimensions uses the argument principle for 
a, see Corollary I4.9l in the present paper. We have not been able to generalize this 
proof either to higher dimensions, even though we think that such a proof may 
not be completely out of reach. At least we have computed the relevant mapping 
degree of a in higher dimensions, and this result, Theorem 14.81 may be considered 
to be the main result in the second part of the paper. Another avenue which we 
have pursued to some extent is the investigation of the general behaviour of the 
curve 7. This curve starts out at the origin and reaches (<957) + only at points where 
the largest inner ball centered at the origin touches (dil) + , see Proposition 14.101 
If we could establish for example that 7 were a smooth curve (no branchings) a 
proof of the inner ball property would not be far away. 

Thus we hope that the partial results we obtain in this paper will turn out to 
be useful in forthcoming attempts to prove the inner ball property for quadrature 
domains, explicitly formulated in Conjecture 14. f I 
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2. Differential geometric preliminaries 
2.1. Notations 

In this section we decompose the coordinates in R™ typically as (u, v) where u = 
(ui,...,u„_i) £ R n_1 , u e R. We denote by (R")+ = G R™ : w > 0} the 

upper half space. Let g be a positive, twice continuously diffcrcntiable function 
defined in a bounded domain D C R™ -1 and assume that g(u) — > as u — > dD 
(the boundary as a subset of R" -1 ). It is convenient to set g(u) = for u ^ D, 
and then g is defined and continuous on all R n_1 . Set 

fl = {{u, v) e R™ : u £ L>, \v\ < g(u)}, (2.1) 

and for any subset A C R™, A+ = A (1 (R™)+. 

The normal line of (c?0) + at a point (u,g(u)) is given in parametrized form 

as 

t ^ (u,0(u))+t(Vs(u),-l)) 
and it intersects the hyperplane R" _1 for t = g(u), i.e., at the /ooi pomi 

p(u) =u + g(u)Vg(u). 

Let 

7V„ = {(u, (/(«)) + i(Vff(u), -1)) : < i < </(«)} 
be the part of the normal line which is between the base point on (<90) + and the 
foot point. The length of N u is 



\N u \=g(u)^l + \Vg{u)\\ 

For any point x £ + there is at least one closest point (u,g(u)) on (dil) + . 
Then x S N u and, in particular, x £ B(p(u), \N U \). Note that N u is one of the 
radii in B(p(u), \N U \). It follows that we always have inclusions 

n+ c |J N u c |J B(p(u), \N U \). 
ueD ueD 

Also, 

Oc |J B(p(u),\N u \). (2.2) 
ueD 

2.2. The first and second fundamental forms 

We shall discuss (dil) + from a differential geometric point of view. We consider it 
as a Riemannian manifold of dimension n — 1 and with coordinates u\, ...,u n -\. 
The Riemannian metric is that inherited from R™. We shall write certain quantities 
considered as vectors in R" in bold. These involve the "moving point" on (<90) + 
(or lift map D (90)+) 

x = x(m) = (u,g(u)), 
its differential (a vector- valued differential form) 

dx = (dui, du n -i, ^ — dui), 

i=i ° Ul 
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the unit normal vector 

(Vg(u),-l) 
n = n (u) = — , 

VIVgHI 2 + i 

and its differential dn (which becomes somewhat complicated when written out in 
components). For a point x G + we sometimes write it in bold if we think of it 
as the vector from the origin to x. 

A tangent vector £ on (cMl) + may be thought of in an abstract way as a 
derivation: 

n ~l o 

4 — 1 

Letting this £ act on the moving point x(u) or, equivalently, letting the differential 
dx act on £ gives the same tangent vector regarded as a vector embedded in M n : 

n—l ^ 

(dM) = (6,...,&-i,E&a£)- 

i—l 1 

In classical differential geometry (see for example [3]) one associates to any 
hypersurface in Euclidean space two fundamental forms. The first fundamental 
form is the metric tensor, which gives the inner product on each tangent space. It 
is in our case 

ds 2 — dx • dx = du\ + ■■■ + du 2 l + (S~* —^-duA 2 

dui 

i—l 

E, r og og , 

where the dot denotes the scalar product in M™. The second fundamental form is 
(up to a sign) 

1 d 2 g 
dx ■ dn — — = > — — - — dui ® du^. 

It contains information on how n rotates in different directions, i.e. on how (dVl) + 
is curved within R". The above expression for dx • dn can be derived as follows: 
Since n is orthogonal to (dfl) + we have J^- • n = for all i. Therefore 

dx dn d , dx . d 2 x „ ,„ d 2 q , 

-(— • n) - — — — • n = - (0, ■ " ■• 



dui duj duj dui duiduj duiduj 



1 d 2 g 



y/1 + |V.g| 2 dmduj ' 

proving the formula. 

Both fundamental forms are symmetric covariant 2-tensors, i.e., symmetric 
bilinear forms on each tangent space. When considered as bilinear forms we shall 
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call them A and B respectively (namely A = cfac • dx, B = dx ■ dn). The above 
formulas then mean that for (abstract) tangent vectors £ and r\ we have 

n— 1 n— 1 r\ o 

i—l 3,j=l 

Note that A is positive definite. 

The principal curvatures of (dil) + and the corresponding principal directions 
are the eigenvalues and eigendirections of the derivative of the Gauss normal map, 
taking any point x 6 (dft) + onto the normal vector at the point: n(x) 6 S"™ -1 ; 
the derivative is the induced linear map between the tangent spaces, intuitively 
dx i->- dn. In terms of abstract tangent vectors this becomes C : £ i-> rj, where 
(dx, rj) = (dn, £}, and one can also say that it is the map obtained when expressing 
the second fundamental form in terms of the first: 

A(C,C£)=B(C ) t). 

It is well-known (and easy to prove) that these eigenvalues and eigendirections 
coincide with the stationary points and corresponding stationary directions for the 
quotient of the fundamental forms, namely for the map 

(£ G R n_1 , £ 7^ 0). For example, the smallest eigenvalue (the smallest principal 
curvature) coincides with the minimum value of B/A. 

2.3. The function $ 

Next we introduce the function 

Mu) = \(\uf+ g (uf)= l -\x(u)\\ (u€D) 

namely half of the squared distance from points on (<9£1) + to the origin. More 
generally, for any c = (a, b) with a € R™^ 1 , b > we set 

*c(«) = \{\u a| 2 + (g(u) bf) = l|x(u) - c\\ 

considered to be defined for those u € D for which g{u) > b. Then $ c is twice 
continuously differentiable with 

— = Ui -a i + (g(u)-b) — , 
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<9 2 $ c dg dg d 2 g 

= % + —— + (9(u) - b) (2.3) 



dmduj dui duj duiduj 

<9 2 $ d 2 q 

b^~- (2-4) 



duiduj duiduj 



For c = we can write 



duiduj dm duj + |Vg| 2 duiduj 

Thus considering the Hessian matrix V 2 4> = ( g ®.g u . ) as a tensor, or bilinear form, 
we see that it is related to the two fundamental forms by 

V 2 $ = rfx • dx + \N u \dx ■ dn = A + \N U \B. (2.5) 

We finally notice that $ is related to the foot point map p by 

V$ = p. (2.6) 

More generally, for any a G ffi.™ -1 we have 

V#(„,o)(«)=p(u)-o. (2.7) 

2.4. The Poincare metric 

The Poincare metric in (R") + is given by 

n— 1 

The geodesies with respect to the Poincare metric are the vertical straight lines 

{u = constant, 
v > 

together with all vertical semicircles with centers on R" _1 , namely all curves of 
the form 

I u G L, v > 0, 



\\u- a\ 2 + v 2 = r 2 , 

where L is a straight line in a£L and r > 0. 

We consider now the variable transformation 

T:(u,v)^(s,t), 

where s G K n_1 , i G K, defined by 

s = u, 



(2.8) 



t = I(M 2 +t, 2 ) 



and, conversely, 

I u 



v=y/2t- \s 



2 



(2.9) 
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Thus T is one-to-one and takes (R™) + in the (it, w)-space onto the epiparabola 

W = {( Sl t)eR n :t>±\s\ 2 } 

in the (s,<)-space. By the definition of $, T maps fl + onto the set 

V = {(s,t) : s e D,t < $(«)}, 

and it maps the graph of g onto the graph of <E>. If <I> is extended to all R™ -1 by 
setting $(it) = \\u\ 2 outside D (this corresponds to extending g by zero outside 
D) then W \ V is the epigraph of $. 

Lemma 2.1. The above map 

T : (R™)+ -y W 

gives a one-to-one correspondence between the set of geodesies with respect to the 
Poincare metric in (R") + and the set of straight lines in W . 

Proof. Let 7 be a geodesic in (K™) + . If 7 is a vertical line, then also T(~f) is 
a vertical line in W, hence a straight line. If 7 is of the form (|2.8|) . then just 
performing the substitution (I2.9[) we see that T(^f) becomes 

[t = s -a+ i(r 2 - \a\ 2 ), 

hence it is a straight line. And the arguments can easily be run in the other 
direction: every straight line in the (s, i)-space which has a nonempty intersection 
with W is of the form T(j) for some geodesic 7 in (R™) + . □ 

In the new coordinates (si, . . . , s«-i, t) the Poincare metric takes the form 

1 n— 1 n—1 

ds2 = f2f-l,Rg ( (2t - |S|2) ^ ( ^ )2 + ( ^ S ^ Si)2_ 

^ 1 1 ' 3 = 1 3=1 

n-1 

- Sj {dsj (g>dt + dt(g> dsj) + (dt) 2 ) , 
3=1 

as an easy calculation shows. Hence this is a metric in W which is complete and 
has constant negative curvature, and whose geodesies are Euclidean straight lines. 

3. Equivalent criteria for the inner ball condition 

Below we state equivalent criteria for the domain defined by (12.11) to satisfy the 
inner ball condition for balls centred on the symmetry plane. The definition of the 
inner ball condition can be taken to be statement (v) in the theorem. 

Theorem 3.1. With assumptions and notations as in Section^ the following state- 
ments are equivalent. 
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(i) The restriction of $ to any convex subdomain of D is convex. (Note that D 
is not required to be convex itself.) Equivalent formulations: the matrix of 
second derivatives is positive semidefinite: V 2 $ > in D; the extension o/$ 
to M 71 is convex; W \ V is convex as a set. 

(ii) For every c — (a,b) with a £ R" -1 , b > 0, the function $ c is convex (alter- 
natively: strictly convex) when restricted to any convex subdomain of the set 
of u € D for which g(u) > b. 

(iii) The restriction of the foot point map p to any convex subdomain of D is 
monotone, i.e. 

(p(u)-p(u')) -(u-u')>0 
for all u, u' in the subdomain. (Here the dot denotes the scalar product in 
M™ _1 .) Equivalently, the matrix (^- + jf^f-) is positive semidefinite. 

(iv) 

= |J B(p(u),\N u \). 
ueD 

(v) There exist radii r — r(u) > such that 

n=\J B(u,r(u)). 

ueD 

(vi) 

N u R N u > = for 11/11' (u, u e D). 

(vii) The principal curvatures at any point (u,g(u)) of (dft) + are > — r^-r- 

(viii) Every point in Q + has a unique closest neighbor on (dfl) + . 

(ix) Every point on (dQ) + is a closest point on dfl for some point in D. 

(x) \ Q is convex with respect to the Poincare metric in 



Proof. (i)=Kii): Fix 6 > and ueD. For any £ E E n_1 with |£| = 1, let 

d 2 g 



" _1 d 2 g 
^ du.dufi^' 

i,j=l 



If a > — „, A — rv then it follows from (12.31) that 

— 2(g(u)-b) 



9 2 $r ^ 1 



.4^, duidu 

1,3=1 J 

while if a < — 2 ( g (u)~b) equation (|2.4[) shows that 



f-r± x duiduj 2(g(u) - b) ' 

From these two inequalities the strict convexity of $ c follows. (We even get a 
uniform lower bound of V 2 $ c in {u £ D : g{u) > b}.) 

(ii)^(i): Just let b -> in 
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(i)o-(iii): This follows from (|2.6p . which also shows that the matrix (g^-) is sym- 
metric itself (p- = t£$-). 

v OUj OUiOUj ' 

(i)^(vii): As remarked in Subsection 12 . 2 1 the smallest principal curvature of 
coincides with the minimum value of the quotient B /A between the two fundamen- 
tal forms. Therefore it is immediate from (|2.5I) that (vii) is equivalent to V 2< 1> > 0, 
i.e., to (i). 

(i)=> (iv): By (|2.2p we only need to show that 

B(p(u),\N u \) eft (3.1) 

for all u G D. 

Fix any uq G D. Then taking a — p{uo) in (|2.7|) we see that the map u i— > 
^ > (p(uo),o)( u ) nas a stationary point at u — uq. In view of the interpretation of 
^(p(uo),o) as half of the squared distance from (p(ito),0) to (u, g{u)) it follows that 
this stationary point is a (global) minimum if and only if (13. ip holds for u — uq. But 
when $ (equivalently Q( p (u ),o)) is convex then every stationary point of &( p (u ),0) 
is a global minimum. Hence (13.11) holds for u — uq. 

(iv) =^(v): If the representation in (iv) holds then we get a representation as in (v) 
by adding small balls B(u,r(u)) C f2 for those u G D which are not in the range 
of p. (Clearly p maps D into D when (iv) holds, but it need not be onto.) 

(v) =^(iv): Let (u,g(u)) G (9f2) + . Then, if (v) holds, there exist points dj G D and 
Xj G B(a.j,r(aj)) such that Xj — > (u, g(u)) as j —> oo. The smoothness of (<9f2) + 
and the inclusions B(a,j,r(aj)) C f2 imply that aj — > p(u) and r(a,j) — > |iV u | and 
we conclude that B(p(u), \N U \) C f2. Now (iv) follows from (|2.2[) . 

(iv)<=>(ix): Note that (u,g(u)) G (cT2) + is a closest point of a G D if and only if 
a = p(u) and B{p(u), \N U \) C f2. Thus a is determined by (u,g(u)) and it follows 
immediately (in view also of (|2.2[1 ) that all points on (dfl) + are such closest points 
if and only if the inner ball condition in (iv) holds. 

(iv)=^(vi): Assume that (vi) fails, so that there exists a point x G N Ul n N U2 
for some u\,U2 G D, Ui ^ U2- Without loss of generality \{u\, g{u\)) — x\ < 
(u 2 ,g(w 2 )) —x\. Then (ui,g(ui)) G B(x, \(u 2 ,g(u2)) — x\). Since x G N U2 we have 
B(x, \(u 2 ,g(u 2 )) - x\) C B(p(u 2 ), \N U2 \)L){{u 2 , g(u 2 ))}. It follows that (ux,g(ui)) G 
B(p(u 2 ), |iV U2 |). But (ui,g(ui)) f fi, hence (iv) does not hold. 

(vi) =>(viii): If x G fi + has two closest neighbours (ui,g(ui)),(u 2 ,g(u 2 )) G (d£l) + 
then x G iV Ul n N U2 . 

(viii)=>(ii): This conclusion is somewhat related to [S], Theorem 2.1.30 (attributed 
to Motzkin), but we give an independent proof. 

Assume that (ii) fails and we shall produce a point c G fl + with at least two 
closest neighbors on (dfl) + . We may assume that $(o,&) is not convex for some 
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b > in a convex subdomain of Db = {u G D : g(u) > b}. We extend $(o,b) to all 
M™" 1 by setting 

<|, - J $ (o,fc) for D b 



4|m| 2 £br«GR"- 1 \X>6 

1„ ,2 



2V H j +i(fl(«)-6)+n. 

In the latter expression the plus subscript denotes positive part and g is assumed 
to be extended by zero outside D. 

The function < &(o,fc) is easily seen to be continuously differentiable and by 
assumption it is not convex. Let VP be the convex envelope of $(0,6) > which can be 
defined as 

= sup{L(u) : L is affine and < $(o,&)}> (3-2) 
Equivalently, \& is the function whose epigraph, epi* = {(u,v) el" : v > &(u)}, 
is the closed convex hull of the epigraph of 3>(o,6)- I n the present case, because of 
the strict convexity of $(o,&) outside a compact set, the convex hull will be closed 
right away (before taking the closure). Thus 

epi* = conv(epi$ (0jb) ), (3.3) 

conv denoting "convex hull" . 

Since $(o,&) is not convex there exists uq € R 11 ^ 1 such that 

*(«o) < i (0(6 )(«o). 

It is easy to see that the supremum in (I3.2p is attained for each fixed u. Hence 

*(«o) = L(ua) 

for some affine L satisfying 

L(u) < $(o,6) («) for all u. (3.4) 

By p.3| the point (uo,L(uq)) = (uo,^(uo)) S epi\l/ can be written as a 
finite convex combination of points in the epigraph of $(o,&)' Since this convex 
combination must take place within the graph of L, uq is a convex combination 
of finitely many points u (at least two are needed) for which L(u) = $(0^) (it). It 
follows that equality is attained in (|3.4[) for at least two different u. 

Now write 

L(u) = a ■ u + k, 

where a G K n_1 , k G M and the dot denotes scalar product in M™ -1 . Then (I3.4[) 
can be written 

|u - a| 2 + \(g(u) - &)+| 2 > 2k + \a\ 2 , (3.5) 
where equality is attained for at least two different u G K™ _1 . 

Since the left member in (|3.5j) is > with equality in at most one point 
(namely u = a) it follows that 

2fc+|a| 2 >0. (3.6) 
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This positive minimum value can not be attained when g(u) < b because for those 
u the derivative of the left member in (I3.5[) is 2 (u — a) ^ 0. Note that u = a does 
not satisfy (|3.5[) when g(u) < b because of (|3.6[) . 

Thus g(u) > b whenever equality holds in (|3.5p . Clearly also g(a) > b by (|3.5[) 
and (|3.6[) . showing that a G Dj,. It now follows that 



for all u G D with equality for at least two different u satisfying g{u) > b. 

In conclusion we have produced a point c = (a, 6) £ il + having at least two 
closest neighbors on (dil) + . 

(i)'O-(x): That (R™) + \ £7 is convex with respect to the Poincare metric means by 
definition that every Poincare geodesic 7 in (R™) + meets ft at most once (i.e., 7\fi 
has at most one component). By Lemma |2. II this occurs if and only if the straight 
line T(j) in W meets V at most once. This is the same as saying that W \ V is 
convex as a set, which is easily seen to be equivalent to the convexity of $ (W\V 
is the epigraph of the extension of $ obtained by setting g = outside D). This 
proves (i)<^>(x). 



4. Boundaries of quadrature domains 
4.1. Notations 

In the rest of the paper we shall write points in W 1 = M n_1 x R typically as 
x = (x 1 , x n ), with x' = (xx, . . . , £ rl -i) S R" _1 . The letter u will be used to denote 
a specific function, see below. 

Let fi be a finite positive (and not identically zero) Borel measure with com- 
pact support in the hyperplane R™ _1 = {x £ R" : x n — 0}. Then there exists a 
unique quadrature open set f2 C R n for subharmonic functions with respect to /1, 
see [10], [6]. Thus n(R n \ 0) = and 



for all integrable (with respect to Lebesgue measure to) subharmonic functions h 
in fi. By uniqueness is symmetric with respect to the hyperplane R™ -1 . We shall 
assume that ft is connected (i.e., a quadrature domain) and discuss the curvature 
of (9f2) + . In the nonconnected case the discussions will simply apply to each 
component of fl. 



\u - a\ 2 + \g(u) - b\ 2 > 2k + \a 



□ 




Let 



u = U"-U [ 



On the curvature of some free boundaries in higher dimensions 



13 



be the difference between the Newtonian potentials of /i and Q, the latter consid- 
ered as a body of density one. Then 

'u>0 inR", 

< u = outside (4.1) 
^Au = xn-H inR". 

From gU it follows that u £ C 1 ((R")+), hence Vw = on (R")+ \ n (since u 
attains its minimum value there). 

For proving the inner ball property it is enough to consider the case that /z 
has a continuous density on R ra_1 : 

dfi = fdxi . . . dx n -i, (4-2) 

with / = f(x') > on D = fl n R" _ \ / = on R™" 1 \ D and / continuous on 
K n_1 (see Subsection 14.51 for a motivation of this). In that case, u is continuous 
across R™ -1 while -J^- makes a jump of size / across R' i_1 (this is the meaning 
of the distributional Laplacian in (|4. 1[) V In summary, assuming (14.21) u satisfies in 
f2 + the overdetermined system 

Au=l in Q+, 

u=\Vu\=0 on (dn)+, (4.3) 

where in the last equation should be understood as boundary values from the 
upper half space. If \x is not of the form (|4.2p . the system (|4.3p still holds with the 
last equation appropriately reformulated. 

Returning to the general case, the components of Vu are harmonic functions 
in + . By applying the maximum principle to it follows that 

u > in + , 
^— < in n+. (4.4) 

OX n 

Thus u is strictly decreasing with respect to x n , showing that (dft) + is the graph 
of a function, i.e., is of the form x n — g(x') for some positive function g defined 
in D = fl nR n_1 . With a more detailed analysis one can show [4], [I], [6] that the 
function g is actually real analytic and tends to zero at dD. It follows that we are 
in the setting of Sections [2] and [H 

The aim of the forthcoming analysis is to approach the following conjecture. 

Conjecture 4.1. Any subharmonic quadrature domain as above satisfies the 
equivalent conditions in Theorem \3.1\ 

As have already been remarked, Conjecture 14.11 has been settled in the case 
n = 2. Two completely different proofs appear in [7] and j8]. A third proof, per- 
haps somewhat related to the one in [8] , will be given in the present paper (Corol- 
lary H~<iJ). 
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4.2. The function p and the 2-form uj 

Let (r, 9) be spherical coordinates in R™ in the sense that r = |x| = x\ + ■ ■ ■ + x^, 
6 = ^ G (if x ^ 0). Then r-jjp = x • V = n^- + • • • + r B ^. In order to 

analyse the curvature of (<9f2) + we introduce a function p and a 2-form u> according 
to 



1 o du , N , 
p = -r 2 - r— - (n - 2)u, (4.5) 

CJ = r<ir Adu= d(-r 2 du). 



Thus, 



lj = uiijdxi A efej = ujijdxi ® dxj 

(cfei A da;^ = dajj (g> dxj — dxj ® ) with 

du du 

(*>ij = Xi- xj — . (4.6) 

dxi 

Note that the restriction of p to (<9f2) + coincides with the function $ in Subsec- 
tion [53] when the latter is regarded as a function of x' = (a?i, . . . , x n -i)'- with 
(dtt) + given by x n — g(x') we have 

$(x')=p(x',g(x')). 

Recall that the Hodge star operator [3] is defined on basic forms 



as 



dxi ± A • • • A dxi 



*a = ±dxi p+1 A • • • A dxi n , 



where the sign is chosen so that 

a A *a — dx\ A • • • A dx n = dx (the volume form) 

and (ii, . . . , i n ) denotes any permutation of (1, ... , n). Interior multiplication of a 
p-form a with a vector field £ is the (p — l)-form obtained by letting £ occupy the 
first position when a is viewed as an antisymmetric operator acting on p vector 
fields: 

(i(€)«)(C 2 ,...,C P ) = a«,€ 2 ,...,0- 

We shall also need the Lie derivative with respect to a vector field £, and recall 
that its action on any p-form a is 

L $ (a) = i(£)da + d(i{€)a), 

and that its action on a vector field rj is the commutator 

L^r] = [£,*/], 

when all vector fields are considered as differential operators (£ corresponds to the 
directional derivative £ ■ V etc.). 



On the curvature of some free boundaries in higher dimensions 



15 



In our Euclidean setting of differential geometry there are natural one-to- 
one correspondences between vector fields, 1-forms and (n — l)-forms. The basis 
elements correspond to each other according to 

&i <H> dxi <-> *d,Xi, 

and the forms/vectors simply keep their coefficients with respect to the above 
bases under the identifications. This means for example that if a, b are vectors, a 
a 1-form and /3 an (n — l)-form, then 

a <H» a if and only if i(a)dx = *a, 

bo/? if and only if i(h)dx = f3. 

The definition of oj can be expressed in terms of interior multiplication and the 
Hodge star operator as 

*ui = —i(r—)(*du). 
or 

4.3. The Cauchy-Riemann system 

The fundamental relationship between p and uj is expressed in the following Cauchy- 
Riemann type system. 

Theorem 4.2. We have 

*dp = d(*uj) in (4.7) 
Written out in components this is 

Z-±% <*-» 



Proof. It is easy to check that (|4.7j) . when spelled out, simply becomes (|4.8|) . so 
we just verify (|4.8|) : Differentiating (|4.6|) we get 



dujkj d 2 u du 
dxj dx 2 dxk 3 dxjdxk 

when j 7^ k. For j = k there is an additional term J^-, which cancels the term 
— J^- . Summing over j therefore gives 

duj k j _ /\u (n ll— — — r——(———*) 

t—' dx-j dxk dr dxk 

j=i j 

On the other hand, differentiating (|4.5[) gives 

dp . . fti d . du 

— = x k -(n- 2)- — (r— ). 

oxk uXk oxk or 

By a change of order of differentiating, and in view of Am = 1 in fi + , (|4. 8[) now 
follows. □ 
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Introducing the coexterior derivative 5, denned on p-forms by 
the system (14.71) can be written 

dp = Slu. (4.9) 

Trivially Sp = 0, and since u> by definition is exact, duo — 0. Recall that the Hodge 
Laplacian is 

A = -(Sd + dS) 

and that it in our Euclidean setting simply equals the ordinary Laplacian acting 
on the coefficients of the forms it is applied to. From (|4.9p we find (using S o S = 
etc.) 

Corollary 4.3. 

(8p = 0, jdu = 0, 

\Sdp = 0, \dSu} = 0, 

in particular 

Ap = 0, Au = 0. 

4.4. The vector field £ 

Next we introduce the vector field 

£ = Vp = x - (n - 2)V« - V(x • Vit) 

d 

= x — (n — 1 + r— )Vu. 

or 

For any function <p, 

i(\Jip)dx = *dtp, 

hence the Cauchy-Riemann system (|4.7|) can be written in terms of £ and ui as 

i(£)dx = d(*ui). 
Immediate consequences, in view of i(£) o i(£) = 0, are 

i(0d(*w) = 0, 
L<:(*u)) = d(i(£)(*u))). (4.10) 

When n — 2, i(£)(*w) = (because it is an (n — 3)-form), hence L^(*u) = 0. 
We proceed with a geometric interpretation of £ on (9f2) + . 

Proposition 4.4. At any point x £ (9J7) + , £ = Vp equals the orthogonal projection 
of the vector x onio the tangent space of (<9f2) + at x. in particular, £ is tangent 
to (dfl)+, i.e., = on (50)+. 
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Proof. Since all components of *lu vanish on (dfl) + , the (n — l)-form i(£)dx — 
d(*ui) vanishes along (<9f2) + (i.e., its integral over any piece of (<95!) + is zero). This 
means exactly that £ is tangent to (dft) + . But £ = x V(x • Vu — (n — 2)u) and 
V(x ■ Vu + (n — 2)u) is orthogonal to (dft) + because x • Vu + (n — 2)u — on 
(dft) + . It follows that £ is the projection of x onto the tangent space of (9f2) + . 

□ 

Corollary 4.5. A point x = (x',x n ) e (<9J7) + is stationary for |r 2 on (dQ) + 
equivalently, x' is stationary for $, if and only if £ — at x. 

Example 1. In two dimensions, with the identification M 2 = C and using (x,y) as 
coordinates (and with z = x + ly) we can define the Schwarz function in Q + by 

oz 

It is holomorphic in f2 + and extends continuously to (dQ) + with 

S(z) = z on 

In addition, by (|4.3p . 

Im S(x + iO) = f{x) on d{Vt + ) n R. 

In the two dimensional case *uj is a 0-form, i.e., a function, and it is exactly the 
conjugate harmonic function of p. More precisely, p, *o> are the real and imaginary 
parts of the analytic function ^zS(z): 

1 r,/ \ 1 2 du ,du . . 
-zS (z = -r -r— + i— = p + i * w. 4.11 

2 2 or ay 

Example 2. In the case of three dimensions we can identify the 2-form ui with the 
vector field 

u = xx Vu = curl (-r 2 Vu). 

v 2 ; 

The Cauchy-Ricmann system *dp = d(*cu) then takes the form 

Vp = curl 

In terms of £ = Vp and u> we thus have the equations 

div u) = 0, 
curl £ = 0, 
£ = curl w. 

In particular div £ = 0. Note the triple curl identity curl curl curl (|r 2 Vu) = 0. 

Notice also that the vector field u) is everywhere tangent to the family of 
spheres |x| = constant. Thus, in spherical coordinates (r, 0, <p), where 9 (now) is 
the angle between the vector x and the X3-axis and <p is the polar angle for the 
projection of x onto the (ari, a;2)-plane, on writing 

u> = uj r e r + uigeg + ui^e^ 



18 



Bjorn Gustafsson and Makoto Sakai 



we have ui r = 0. On D = finM 2 , i.e., for 6 = § we have uj v = -r^ = \rf(x') > 
because / > on D. 

For the interior multiplication, 

= <*>•£ = (x X Vti) • (x - V(x • Vu) - Vu) 
= -(x x Vu) • V(x ■ Vu). 

As to the Lie derivatives of uj (vector), lu (1-form), *uj (2-form) we have, when the 
results are identified with vector fields, 

L € (w) = [e,a;] = (£.V)a>-(wV)£, 

L^(uj) = as a vector = curl (u> x £) = (£ • V)u> — (u> ■ V)£, 
L^(*w) = as a vector = V(£ • u>) = (£ ■ V)u; + (a; • V)£. 

We now return to the general case. Since, by (|4.4I) . Vu ^ in f2 + , the integral 
curves of Vu constitute a smooth (even real analytic) family of curves which fill 
up f2 + . The curves are directed downwards (Vu-e„ < 0), start immediately inside 
(<9£1) + (recall that Vu = on (dfl) + ) and end up at points of suppu. Thus a 
subdomain R C f2 + which is bounded by integral curves of Vu may be thought of 
as a tube going from (dfl) + to D C M n_1 . 

Proposition 4.6. Let R C f2 + be a domain such that dRHfl + is smooth and consists 
of integral curves of Vu. Then, assuming dfi = fdx\ . . . dx n -\ with f continuous 
on D, 

m denoting Lebesgue measure. 

Proof. By assumption, Vu • n = on dR fl f2 + (n the normal unit vector on OR) . 
Therefore, 



m(R) = / Audx = / Vu • nda 
Jr. JdR 



Vu ■ nda + / Vu ■ nda + / Vu • nda 

dRnR™- 1 JdRnn+ JdRn(dU) + 

du \ 
(--—)dxi . . . dx n -i + + = -fi(dR n IET" 1 ). 

□ 

4.5. The curve 7 and the map a 

Next we shall study the set where u vanishes: 

7 = {xeSl + : oj(x) = 0}. 

The vanishing of a; at a point x means by definition that all components of 
ui vanish. Note that, since rdr 7^ and du ^ in f2 + , to vanishes if and only if 
rdr and du are proportional, i.e., if and only if Vu is parallel to the vector from x 
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to the origin, and hence points towards the origin. This is also the same as saying 
that the projection 

Vas r u = Vu — (Vu • e r )e r 

of Vu onto the tangent plane of the sphere dB(0,r) vanishes. Hence 

7 = {i£ fl + : — Vu(x) is parallel to the vector x} (4-12) 

= {ie!l + : S7dB r u( x ) = 0}- 

Recall that x denotes the vector from the origin to the point x. 

A possible proof of Coni ecture 14.11 may be based on establishing that 7 is a 
smooth curve, when nonempty. The strategy would then be as follows. Suppose 
f2 does not satisfy the inner ball condition, i.e., that the equivalent conditions in 
Theorem 13. II do not hold. Then there exists by condition (viii) a point c = (a, b) £ 
fl + which has at least two closest neighbors on the boundary. Here a £ M n_1 , 
b > 0. By "localization" (see [6], [13]) the set !l+={ie!l+:i„> b} will still be 
the upper half of a quadrature domain, namely for a measure with support on the 
hyperplane x n — b. Moreover, this measure will have a continuous density function 
/ on x n = b, namely given by 

f(xi,...,x n -i) = -2^T (xi,...,x n -i,b) 

(cf. (|4.3|) V This is because the localization measure on x n = b can be obtained 
from u by keeping u unchanged in x n > b and defining a new continuation of u to 
x n < b by reflection in x n = b. Note by (|4.4[) that the above / is strictly positive 
at points belonging to f2. 

Now we translate our system of coordinates so that c becomes the origin in 
the new coordinate system, and we also adapt the notations in general to fit the 
new coordinates. In this new situation is a subharmonic quadrature domain for 
a measure [i in the hyperplane x n = with 

d/i = fdxi . . . dx n -x, (4-13) 

/ continuous on R n_1 , />0onD = On R"" 1 , and the largest ball B(0,r ) 
centered at the origin and contained in fl touches (<9f2) + in at least two points. 
The set 7 is expected to meet (dfl) + at these points (cf. Proposition 14. 101 below) . 
and in addition reach the origin, thus it need to branch. Thus if we can prove 
that 7 is a smooth curve we reach a contradiction. The conclusion then is that the 
original fl will have to satisfy the inner ball condition. 

Motivated by the above we assume from now on that fi is of the form (|4.13[) 
with / > on D and continuous on R ra_1 . 

Proposition 4.7. If convfJ, then 7 is empty. If £ then 7 intersects every 
hemisphere (dB(0,r)) + with < r < ro = inf x6 (9n)+ | x |- For r > sufficiently 
small there is exactly one intersection point. At the origin 7 starts out in the 
direction £(0). The closure of 7 reaches f2nR n_1 only at the origin. 
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Proof. To prove the first statement, assume supp/i C {(x\, . . . , S R ra 1 : 

Xi < 0}, for example, and consider the component 

du du 

U\ n = Xl- X n - 

OX n OXl 

of lj. We have Acj ln = in fi+, w ln = on (90)+. On M™" 1 , w ln = 

which by assumption (and (|4.3p ) vanishes for x\ > and is nonnegative (and not 

identically zero) for x\ < 0. Hence u>i n > 0, and so u ^ 0, in f2 + . 

To prove the second statement, consider any hemisphere (dB(0, r)) + , < r < 
ro, and the tangent vector field Vas r u = Vw — (Vw-e r )e r on it. On the boundary, 
<9£>(0, r) nl"" 1 , Vas r M points into the lower hemisphere, because -g^- < 0. Thus 
it follows from well- know index theorems for vector fields (see [3], for example) 
that VdB r u has to vanish somewhere in (dB(0, r)) + . 

Next, assume x,y £ 7 H (dB(0,r)) + (x ^ y). Then x = — AVzt(x), y = 
— /xVu(y) for some A, p, > (x, y denote x and y considered as vectors). An easy 
calculation, using |x| = |y|, shows that 

\Vu(x) - Vu(y)\ = 1 |x-|y| > 1 . nf |x-ay| > l_ 
|x-y| A |x-y| " A «>o |x - y| " 2A' 

By assumption, Vu is smooth in Vt + up to K" -1 and |Vu| > — J^- > c for 
small r > (for some c > 0). It follows that A = |vJ(»| — c> nence A if we let 
r — >• 0. Thus the difference quotients IZhMz^hMI on (<9£?(0,r)) + tend to infinity 
as r — > 0, which contradicts the smoothness of Vu up to M n_1 . The conclusion is 
that there cannot be two different points in 7 n (dB(0, r)) + for r > sufficiently 
small. 

Thus, for small r > there is a unique point x = x(r) £ 7 l~l (dB(0, r)) + . The 
direction as r is (writing x(r) = — A(r)Vu(x(r))) 

x(r) -A(r)(-V M (z(r))) Vu(0) €(0) 

lim — lim = = 

r->o |x(r)| r^o A(r)|V«(x(r))| |Vtt(0)| |£(0)| 

The last statement of the proposition follows directly from the fact that 

"in =xj-§t = -h x if( x ') on P nr_1 ) \ {°}- 

□ 

From the definition (|4.6[) of ojkj we see that 

Ldkj = (XjUJk n — Xk^jn), 

Xn 

hence oj is algebraically determined by the n — 1 components u>i n , ■ ■ ■ ,^ n -i,n- It 
follows that all information of p and u) is contained in the map 

a : tt + ->• M™, .x !->■ (wi„(x), . . . , u) n -i <n (x), p(x)). 

In the case of two dimensions this a = (u>i2, p) is, with a switch of coordinates, the 
same as the function ^zS(z) in (I4.11[) : ^zS(z) = p + i*uj = p + io; 12 . In particular 
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we conclude from the above that 

7 = {x e ft + : wi n (x) = • • • = w„_ 1>n (x) = 0} 
= t7 _1 (the n:th coordinate axis). 
Let J a denote the Jacobian matrix of a: 

Ht ^ J- 1 'I- r: ' ' ' H T 



Ja = 



V 









dx\ 


dx2 


dx n 








dxi 


dx 2 


dx n 


dWn.— -L } n 




duj n - 




dx2 


dx 


Qp 


dp 


dp 




8x2 


dx n 


says that 







J 



tr J a = 0. 

Equivalently, if a considered as a vector field, 

div cr = 0. 

Note also, by Corollary 14.31 that a is a harmonic map (each component of a is a 
harmonic function). 

Clearly a extends continuously to Q + . On (<9i7) + , 

a=(0,...,0,ir 2 ), (4.14) 

hence (<9f2) + , of dimension n— 1, is mapped by a onto a set of dimension at most 
one. It follows that rank J a < 2 at all points of (dtt) + (since a loses at least n — 2 
dimensions). To be more precise, for any vector 77, 



V ■ Vu 2 „ 



77 • Vu„-i,. 



/ 



\ T7 - Vp 

and writing, at a given point on (917) + , 77 as 

*7 = Vt + Vn n 

with r/ t tangent to (<9£1) + and rj n 6 K gives (since w jn = on (9il) + ), 
^°-*7 = fat ' e « + ?7n ^ (n) = (rj t ■ x) e„ + rj n J a (n), 
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where e„ is the n:th unit vector (for the last member Proposition 14.41 was used). 
Thus the range of J a is contained in the (at most) two dimensional vector space 
spanned by e n and J CT (n). 

At points on (d£l) + where \r 2 is stationary, i.e., £ = 0, we have rj t ■ Vp = 0, 
hence rank J a < 1 . Since the trace of J a is zero it follows that all eigenvalues of Jo- 
are zero. In two dimensions it even follows that J a = 0, because J a is symmetric 
when n = 2. 

Returning to the mapping properties of er, we see from (|4. 14[) that a maps 
(dft) + onto the segment \r^ < a n < \r\ of the cr n -axis, possibly except for 
endpoints, where 

ro = inf r, T\ — sup r. 

On D = n n E n_1 we have 

a = {X1 —. ..,x n ^—, -r -r--{n-2)u) 

= {—fix'), -^f{x'), \r 2 - - (n - 2)u). (4.15) 

Thus, since we have assumed that / > on D, a(D) meets the <r ra -axis only if 
G fl. The meeting point is cr(0) = (0, . . . , 0, -(n - 2)u(0)), hence has er„(0) < 0. 

In summary, a maps the boundary d(ft + ) of ft + as follows: £ d(il + ) is 
mapped onto a point (0, . . . , 0, a n ) with a n =-(n- 2)«(0) < 0, 9(fi+) n (M™)+ = 
(9J7) + is mapped to points on the cr„-axis, with a n > > 0, and for the 
remaining points x € D \ {0} of d(il + ), cr(x) does not meet the <7 n -axis. More 
precisely, x and cr(x) are by (|4.15l) located on opposite sides of the cr„-axis. It 
follows from all this that <9(0 + ) encloses the interval [—(n — 2)u(0), |rg] on the 
cr„-axis (— l)"^ 1 times in the sense of degree theory. Thus 

Theorem 4.8. Assume 6 and that fj, is of the form (|4.13[) with f > on D. 
Then the map a has mapping degree (Brouwer degree) {— 1)" _1 with respect to 

each point (0, . . . , 0, t) with — (n — 2)u(0) < t < 

2 

A conclusion is that each value (0, . . . , 0, t), < t < is attained at least 
once, but it can also be attained several times with signatures (signs of the Jacobi 
determinant) adding up to (— l)" -1 . 

In the case of two dimensions a stronger conclusion is possible because the 
Jacobi determinant has constant sign: det J a < 0. Switching the order between u 
and p one gets the analytic function ^zS(z) = p + i*uj (which we still denote o~(z)), 
with nonnegative Jacobi determinant, and all arguments become more familiar, 
e.g., the mapping degree can be identified with the argument variation divided by 
2ir. We repeat everything in this case, with stronger assertions. 

Corollary 4.9. With assumptions as in the theorem we have, in the case n = 2. 
that if € f2 then there is exactly one point on (9£!) + of shortest distance to the 
origin and 7 is an analytic smooth arc from origin to that point. In particular, the 
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equivalent conditions in Theorem \3.1\ hold. In addition, 7 is an integral curve of 

£■ 

Proof. Let ro > be the distance from the origin to (dfl) + . We apply the argument 
principle to the analytic function cr(z) = ^zS(z) = p + i * to in ft + . We have 
er(0) = 0, Imer(V) > for x < 0, Im<r(a;) < for x > 0. For z G (90)+, 
a\z) = ±|z| 2 , hence <r((0fi) + ) C [§r§,oo). 

It follows that as z runs through d(fl + ), then cr(z) encircles each point on 
the open interval (0, \r^) C R exactly once. In addition, no other point on R is 
encircled, although some points lie on the image curve a(d(fl + )). But now 7 is 
exactly the inverse image of R under a, hence it follows that a is univalent on 7. 
This means that t n> er -1 (t), for < t < ^rg, is a bijective analytic parametrization 
of 7. 

Since (<9f2) + is known to be analytic without singular points in the present 
geometry [TU], [T2], [B] the Schwarz function S(z), and hence a(z), extends to be 
analytic in a full neighborhood of (<9il) + . In particular, a is an open mapping in 
such a neighborhood. 

Let zq 6 (d£l) + be a point on distance ro from the origin. Then, with the 
above extended a, ct(zq) — \r\, and any neighborhood U of zo in (M 2 ) + is mapped 
onto a full neighborhood cr(U) of ^r 2 ,. The latter neighborhood contains the in- 
terval {\tq — e, |rg) on the real axis for some e > 0, hence E7 contains the final 
part of 7. Since this is true for any neighborhood U of zq it follows that j!l (dfl) + 
consists of just the point zq. In particular, there is only one point on (dfl) + on 
distance tq from the origin. 

The last statement of the corollary follows from the fact that = (see 

(|4.10p ) in two dimensions, or simply from the analyticity of ^zS(z) = p + i * uj: 
7 is a level line of *u> and £ = Vp is perpendicular to level lines of p, which are 
orthogonal to those of □ 

Now we return to n > 2 dimensions. We still expect that there is only one 
point on (d£l) + of shortest distance to the origin, that 7 reaches (<9£1) + only at 
that point, and that there are no other stationary points of \r 2 . The following 
proposition gives some partial results in this direction. 

Proposition 4.10. The closure 0/7 in R ra intersects (dft) + only at points of (dft) + 
where £ — 0, i.e., where ^r 2 is stationary. If there is a strict local maximum of 
\r 2 at some point x € (dfl) + , then x is in the closure ofj. 

Proof. It is easy to see that, at points in f2 + close to (Sf2) + , V« is essentially 
directed in the inward normal direction. At any point x € (d£l) + where £ ^ 0, the 
vector x from the origin to x has a nonzero angle to the normal vector. It follows 
from (|4.12p that Vu cannot be parallel to x close to such points, hence that 7 does 
not reach points of (dil) + where £ 7^ 0. 
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Next, let x 6 (<9£1) + be a point at which \r 2 attains a strict local maximum. 
Then for some neighborhood N of x 

(TF\{x})nN C B(0,r)DN 

(r = \x\). For a slightly smaller r, say r' < r, S = (0 + n dB(0, r')) n N is a piece 
of a sphere cut off by (dfl) + . In this piece u > 0, while on the boundary u = 0. 
It follows that there is a point on S where u\s attains a local maximum. But at 
such a point the gradient Vtt is perpendicular to S, which means that the point 
belongs to 7. Since r' was arbitrarily close to r, this means that we can produce 
points on 7 arbitrarily close to x. □ 

4.6. The Hessian 

Let 

d 2 u , 



dxidxj 



be the Hessian of u. It is a symmetric matrix with tr H = Au = 1 in The 
rank of H is at least one (since the sum of eigenvalues is different from zero). On 
(<9£1) + the rank is exactly one since Vu = on (dfl) + and hence all derivatives of 
V?i in tangential directions of are zero. 

From r-g^Vu — (x • V)Vu = -ffx, where the last product is matrix multipli- 
cation, it follows that the definition of £ can be written in terms of H as 

£ = x - (n - l)Vu - Hx. 

Thus, on (90)+, 

Hx = x - £, 

and in particular 

#x = x 

at stationary points (for ^r 2 ) on (<9f2) + . 

More generally we can write, for any 77 € M n , 

ii?7 = V(?7 • Vu). 

On (dfl) + 7 r] ■ Vu — 0, hence V(rj • Vu) is perpendicular to (dfi) + . Therefore, on 
(90)+, 

i?t7 = An 

for some A = A (77) € R. It follows that n defines the only eigendirection for H 
corresponding to a nonzero eigenvalue, and since the sum of eigenvalues equals 
one, that 

Hn = n. 

Moreover, it follows that H annihilates every vector tangent to (9f2) + . Thus, on 
(<9S1) + , H simply is the orthogonal projection onto the normal of (<9£1) + : 

Hr] = (rj ■ n)n. 
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Proposition 4.11. With dx = dx\ . . . dx n = dm, 

32 



IQ+ ud -n 

Proof. For any entry (i,j) ^ (n,n), say i = 1, j arbitrary, we have 

/ A" dx i ■ ■ ■ dx » = / d (-K— dx 2 ■ ■ ■ dx n ) 
J n + dxxdxj J n + dx 3 

= f 7T- dx, 2 ■ ■ ■ dx n = 

(note that dx n =0onR n_1 ), while for i = j = n, using the above 

/ ^r^-dx = f Audx = ^m(fi) = ]- [ dfj,. 
Jn+ Cx n j Q+ z z j 



□ 



4.7. Convexity questions 

We finish by some simple observations related to convexity. 

Proposition 4.12. The restriction ofU 11 to any hemisphere {dB{a 1 R)) + with center 
a on W 1 ^ 1 is convex as a function of x' = (x\, . . . , x„_i). 

Proof. This follows from a straight-forward computation of the second derivatives 
of x' i \ U^(x', \J R 2 — \x' — a\ 2 ). Assuming for simplicity that a — and com- 
puting, e.g., the derivatives in the xi-direction in the case n > 3 we have, with 

x n R ~ \ x I i 

My') 



U^x',^R 2 -\x'\ 2 ) = U^x',x n ) 



(n-2)\S n - 1 \ J (|a/_y/|2 + 3-2)^ 

My') 



(n-2)\S n - 1 \J (\y>\2 ~2x> -y* + R 2 )^ 
dU» = 1 f (n- 2) Vl dfM) 

dx x {n-2)\S n - x \} (\y'\ 2 -2x> -y> + R 2 )% ' 
d 2 U^ _ 1 f n(n-2)y 2 dn(y') 



dx 2 (n-2)\S n - 1 \J (\y>\> -2x> -y> + R 2 )^ 
where denotes the spherical area of 5 n_1 . 



>0, 



□ 



Proposition 4.13. Assume that supp ^ C B(a,R) flR" -1 for a ball B(a,R) with 
center a 6 K n_1 . Then the density of classical balayage of \i onto the sphere 
dB(a,R) is convex as a function of x' . 
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Proof. We may assume that a = 0. Set Br = B(0,R). Classical balayage of a 
point mass at any point y e Br onto 8Br gives the mass distribution on 8Br 
with density equal to the Poisson kernel (or normal derivative of Green function) , 
i.e., the function of x e 8Br 

nA ,y> \S n - l \R \x-y\ n \S n - l \R {\y\> ~ 2x ■ y + ' 

It follows that the density /3o of classical balayage of /i is obtained by integrating 
the above with respect to y' G R™ -1 : 

i [ R 2 - \yf , , i\ 

IS^IRJ (\yi\ 2 -2x> -y' + R 2 )% 

Note that the last expression is a function only of x' . One sees immediately 
that as such a function the second order derivative in any direction is nonnegative 
(compare previous proof). In fact, any derivative of even order is nonnegative (a 
similar remark applies to the previous proposition). 

□ 
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